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I .   INTRODUCTION 

The  problem  of  the  reflection  of  an  acoustic  plane  wave  from 
a  plane  interface  separating  two  isotropic  media  has  been  solved 
by  many  physicists  and  geologists.   Exact  forms  of  the  generated 
waves  are  given,  related  to  the  incidence  angle  and  the  acoustic 
impedance$.of  the  media.   Brekhovskikh  extended  the  work  to  deal 
with  a  layer  having  plane  interfaces.   The  nature  of  wave  scatter- 
ed by  a  layer  with  rough  interface,  however,  is  generally  unknown, 

The  backscattering  of  acoustic  wave  from  a  rough  layer  is 
analyzed  in  this  report.   The  model  of  the  layer  considered  con- 
tains a  smooth  interface  in  front  and  a  random  rough  interface 
in  back.   The  Kirchhoff's  approximation  for  evaluating  the  scat- 
tering field  of  a  rough  surface  (Beckmann  1963)  is  extended  to 
deal  with  such  a  layer;  experimental  work  has  also  been  done  for 
this  model  at  different  incidence  angles  and  frequencies. 

A  tentative  try  is  also  made  on  the  evaluation  of  back- 
scatter  from  a  layer  with  a  very  rough  surface  in  the  front. 
Equations  are  derived  for  this  layer  with  no  experimental 
support. 

Owing  to  the  analogies  between  acoustic  and  electromag- 
netic waves,  the  results  of  this  work  can  be  directly  applied  to 
the  same  problems  in  electromagnetic  waves.   The  radar  cross  sec- 
tion is  obtained  just  by  modifying  the  variance  of  the  scattering 
coefficient  with  a  scaling  factor.   Application  of  the  results 
can  be  found  in  the  survey  of  lunar  surface  and  in  geological 
explorations. 
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II.   EQUATIONS  OF  ACOUSTIC  WAVE  MOTION 
2.1  Stresses,  strains ,  and  elastic  constants 

In  an  ideal  isotropic  homogeneous  medium,  a  wave  may  propa- 
gate without  any  loss  of  amplitude  due  to  internal  friction.   When 
the  medium  is  defomable  and  undergoes  a  change  in  configuration 
due  to  the  application  of  forces,  the  body  is  said  to  be  strained. 
(Redwood  1960,  Ewing  1957).   It  is  assumed  that  a  point  P  (X ,    Y, 
z)  is  displaced,  and  the  coordinates  of  the  displacement  are  (u, 
v,  w) .   An  adjacent  point  Q(x+6x,y+6y ,z+6z)  is  displaced  by  (u+ 
6u,v+5v,w+6w) .   By  Taylor's  theorem,  neglecting  the  higher  terms 
under  the  assumption  of  small  perturbation, 

u  +  6u  =  u  +  {H  6x  +  |H  6y  +  |H  6z    , 

v  +  6v  =  v  +  |X  6x  +  jX  6y  +  |Z  6Z     ,  (2.1) 

w  +  6w  .  w  +  |*  6x  +  »S  6y  +  ||  62     . 

The  analysis  may  be  simplified  by  writing 

3u  1Z  -  e  3w  _ 

3x  =:  xx    '     3y  "  yy    '     3z  '"  zz    ' 

3x   3y  "  exy  "  yx        '        3z   3x  "  zx  ""  xz 

3y    3z    yz    zy     ' 

—  -  —  ■  co  iH-iZso)  3w    8v  m       | 

3x  "  3y     z    '   3z    3x     y    '    3y  "  3z  "'   x 
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Then  the  displacements  may  be  rewritten  as 

u+6u=u+(exx6x+iexy6y+iexz6z)  +  (u)y6z-u>z5y) 

v+6v=v+(^eyx6x+eyy6y+|eyz6z)  +  (a)26y-u)x6x)     ,         (2.3) 

w+6w=w+(^ezx5x+jez  6y+ezz6z)  +  (ujx6y-(o  6x) 

e   ,  e   ,  e   represent  simple  extensions  of  the  medium  in  the 
xx '  yy   zz   r 

vicinity  of  P  (x,  y,  z) ; 

e   ,  e   ,  e   represent  the  shear  strains;  and  o>  ,  ■  ,  w  repre- 
xy '  yz'   zx  *  x   y   z 

sent  the  rotation  of  the  element  as  a  rigid  body. 

To  express  the  displacement  in  vector  form 

s  =  ux„  +  wz^  +  vy„     ,  (2.4) 

o     o     o 

where  x   .  v   /  and  z  are  unit  vectors  in  the  directions  of 
o    o        o 

x-,  y-,  and  2-   axes,  respectively. 

For  small  perturbation,  the  shear  strain  is  so  small  that  it 
has  no  contribution  to  the  volume  change.  Then  the  cubic  dilata- 
tion A  is  defined  as 

A  =  lim        (6x+exx6x) (6y+eyy6y) (6z+ezz6z)r6x6y6z  {25) 
6x,6y,6z->-0  6x6y6z 


=  e   +  e   +  e 
xx    yy    zz 


=  div  s    , 
and  the  rotational  displacement  u  becomes 


-4- 


•*■«*■»■      ■*• 
x  o    y-^o    z  o 

=  j  curl  s    .  (2.6) 

To  discuss  the  force  acting  on  an  element  of  volume  in  a 
medium,  nine  components  of  stress,  which  have  the  units  of  force/ 
area,  are  required.   Let  the  stresses  be  denoted  by  t *  . ,   where 
the  first  subscript  is  associated  with  the  axis  normal  to  the 
place  on  which  the  stress  acts,  and  the  second  subscript  to  the 
direction  of  the  stress.   In  an  isotropic  homogeneous  medium, 
there  are  only  two  elastic  constants,  Lame's  constants  X  and  y . 
to  relate  the  stress  and  strain  together.   In  the  case  of  ideal 
fluid  and  ideal  gas  u  ■  0.   The  stress  and  strain  relations  by 
Lame  are  as  follows: 


t  =X  (e  +e  +e   )+Zye  =XA+2ue 
xx   v  xx  yy  zz'  ^  xx       xx 


t  «x  (e  +e  +e   )+2ue  =XA+2ye 
yy    xx  yy  zz     yy     H  yy 


/ 


t  =X  (e  +e  +e  )->2ue  =XA+£ue      , 
zz    xx  yy  zz7   M  zz    fcH  zz    ' 

t   =t   =pe  =ye      ,  (2.7) 

xy  yx   xy   yx    ' 

t  =t  =ue  =ue      , 
yz   zy  M  yz  M  zy     ' 


t  =t  =ue  =ue 
zx  xz  M  zx   xz 


From  (2.7),  other  elastic  constants  can  be  derived  as  follows, 
(Lamb  1925)  : 
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(a) .   Uniform  stress  and  dilatation,  t v  =t   =t   ,  and 


x       =T   -T   =(X+|-u)A=kA 
xx  yy  zz    3 

2 
k=(X+-=-y)  is  the  compressibility, 


(b) .   Shear  stress  t    as  defined  by  Eq.  (2.7) 

xy 

T        T        T 

xy  _  yz  _  zx 

M   e     e     e 

xy    yz    zx 


is  the  coefficient  of  rigidity, 
(c) .   Longitudinal  stress,  t  v=tzz=0 

m    y(3X+2y)      m 
xx      X+y    xx     xx    ' 


E  is  the  Young's  Modulus. 
(d) .   Poisson' s  ratio,  a    ,  represents  the  ratio  of  lateral 
contraction  to  longitudinal  extension.   In  this  case, 

Tyy  =  Tzz  "  °   '   eyy  "  ezz  =  oexx   '  and  °  =  aTI+yT 

The  constants  X,  y,  k  can  be  expressed  in  terms  of  Young's 
Modulus  and  Poisson 's  Radio: 

X  ■  Lame's  constant  = 


V 


(l+o) (l-2a) 


=  Rigidity  =  ^1+p)     ,  (2.8) 


E 
k  =  Compressibility  =  ^  /i_j  \ 


2.2   The  equation5Qf  motion 


According  to  Newton's  equation  of  motion,  the  x-component  of 
the  resultant  force  on  the  volume  element  is 
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(p6x.6y5z)0  =  (txx  ♦  ^x  6x-Txx)6y62 


+  <Tyx  +  -^*6y-Tyx)6Z6x 


3tzx 


+  ^«  +  -arl8M«,6x6y 


32u  .  ,  3xxx  ,  3Tyx  t  3tzx  , 
or'     p  ap-"  (  "ax"  +  "ay"  +  ^i~  > 


/ 


where  p  is  the  density  of  the  medium.   For  an  isotropic  homo- 
geneous medium,  Equation  (2.7)  holds,  giving 


"  0  ■  fc  (u+2ueXx»  +  37  (UV  +  h  <»•«> 


5A         3to     3u> 


Similar  arguments  apply  to  the  y-,  and  z-  components,  so  that 

32s  -*■ 

p  jp-  =  (X+2y)  grad  A  -  2p  curl  w     .  (2.9) 

Equation  (2.9)  is  the  equation  of  motion  in  an  isotropic  medium. 
It  is  usually  rewritten  in  the  form  of  displacement  potentials  <J> 
and  %   by  the  following  transformations 

s  =  -  (grad  <J>  +  curl  $)     ,  (2.10) 

div  ty   -   0 


-7- 


Substituting  (2.10)  into  (2.9), 

p — -(grad  <j)+curli^)  =  (A+2y  )grad(V<()) -ucurl  (curlcurl^)  , 
at 


p— 9(grad<{i+curl|)  =  (X4Z;M)grad(V    <j>)+ucurl(V   J),  (2.11) 

3tz 

2  .  2 

where   V     denotes   the  Laplacian,   defined  by  V  ♦  ■  Div.    Grad   $ 

_  2-*-  •*■  ,    •*■    ,  ■+ 

for  a   scalar   <|> ,    and   V    $   =   Grad  Div   ijj-curl   curl   ij)   for   a  vector   ij/. 

In  the  cartesian  coordinate   system, 
2  2  2 

v2=i-  +  i-  +  i- 

2      2      2 

3x     3y     3z  • 

By  taking  the  divergence  and  curl  on  Equation  (2.11)  the  scalar 
and  vector  potentials  can  be  separated: 


924>   ..   r    2    2  2   .     A  +  2p    .  ,-    ,_. 

— £■  -  C.    7    <f>    ,    C.      =  r    ;  (2.12) 

3t  <  ^  p 

i4-   C    Vj    ,    C   2   =  H.      .  (2.13) 

3t2  s  s  p 

In  Equation  (2.11)  and  (2.12),  the  scalar  displacement 
potential  $  travels  with  velocity  C»  and  involves  no  rotation;  it 
is  called  "longitudinal",  or  "compressional" ,  or  "dilatational" 
or  "irrotational" ,  or  "P-Wave";  the  vector  displacement  potential 
i|»  travels  with  velocity  C   and  involves  rotation;  it  is  called 
"shear",  or  "transverse",  or  "lateral  "  ,  or  "equivoluminal" ,  or 
"rotational",  or  "distortional" ,  or  "S-Wave".   The  names  Longi- 
tudinal and  Shear  waves  shall  be  used  throughout  this  report. 
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In  ideal  gases  and  fluids,  the  ordinary  acoustic  wave  init- 
iates a  wave  motion  in  which  the  sign  of  dilatation  changes  very 
fast  so  that  there  is  no  time  for  sensible  transfer  of  heat  be- 
tween adjacent  portions  of  the  medium.   The  flow  of  heat  hardly 
sets  in  from  one  element  to  another  before  its  direction  is  re- 
versed, and  the  conditions  are  practically  adiabatic.  Moreover/ 
since  u  ■  o,  no  shear  wave  propagates  in  an  ideal  gas  or  fluid. 
The  Lame's  constant  X  for  an  ideal  gas  is  found  to  be 

X  =  K  ,  =  \Pn  (2.14) 

ad     o 

K  ,:  adiabatic  compressibility, 

X   :  (specific  heat  at  constant  volume) /(specific 

heat  at  constant  pressure) , 

P   :  gas  pressure. 

In  an  ideal  fluid,  X  =  K  ,  is  usually  influenced  by  many 

factors.   X  is  related  to  the  longitudinal  wave  velocity  C& 

by  C,  ■  /x/p  ,  and  an  empirical  equation  giving  C,  in  water  is 

C^  =  141, 000+4. 21t-3.7t2+110s+0.018d,      (2.15) 

C.  =  longitudinal  wave  velocity  (cm/sec) , 

t  =  temperature   (C) , 

s  =  salinity   (1/1000  in  weight) , 

d  =  depth   (cm) . 
2.3   General  equation- for  damped  waves 

The  dissipative  forces  in  acoustic  wave  motion  are  propor- 
tional to  the  velocity  of  the  particles  in  the  medium  for  small 
perturbations.   The  force  of  the  elastic  stress  must  both  accel- 
erate the  medium  and  overcome  the  dissipative  forces,  so  the 
equations  of  motion,  Equations  (2.11)  and  (2.12)  are  modified 
(Kinsler  and  Frey  1950)  as  follows: 
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^4  +  7   if-  C,2V24>,  (2.16) 

9t2    p  3t     * 

ll  +  a  ||  =  C  V*  ,  (2.17) 

3t2    p  3t     S 

,.1.        c    (mass)  (time) 
where  n  has  the  unit  of    (volume) 

The  general  solution  of  Equations  (2.16)  and  (2.17)  has  the 
usual  form  E  =  Ae1  *X  •r_U)t)  ,  where  r  is  the  distance  from  certain 
reference  point.   Substituting  into  Equations  (2.16)  or  (2.17), 
the  propagation  constant  X  is  found  to  be 

X  =  K+Kat 
where 

K  =  o)/C   is  the  wave  number, 
a  =  n/2pC  is  the  attenuation  constant. 
For  water  at  room  temperature,  the  attenuation  constant  is 
about  one  thousandth  of  the  value  of  air.   For  this  reason, 
ultrasonics  are  very  well  suited  for  underwater  signal  transmis- 
sion as  opposed  to  the  case  in  air.   As  for  the  possible  shear 
waves  transmission  in  viscous  gas  and  liquid,  the  skin  depth  is 


found  to  be  /2£/po>  ,  where  K  is  the  shear  viscosity. 

These  waves  behave  very  much  like  electromagnetic  waves 
penetrating  a  metal.   For  water  at  1  mc,  the  skin  depth  is  only 
5.6  x  10~  cm.   In  most  cases,  shear  wave  propagation  can  be 
neglected  in  gases  and  liquid. 


3.  REFLECTION  OF  PLANE  WAVE  FROM  A  PLANE  INTERFACE 

3.1  Boundary  conditions 

The  continuity  of  normal  and  tangential  components  of  dis- 
placements and  stresses  across  the  interface  give  the  following 
boundary  conditions  (Brekhovskikh  1960)  : 

-*■        ,-*■-*■   \ 

n  .  (s1~s2)=o  ; 

n  x  (s1~s2)=o  ;  (31) 

n  .  (t1~t2)=o  ; 

n  x  (t,-t2)=o  . 

Taking  the  z-axis  as  the  normal  to  the  plane,  Equation 
(3.1)  can  be  expressed  in  terms  of  displacement  potentials  $ 
and  <(i  as, 

n   s  w      vaz        3x  3y   '    ' 

nxs=-uxo+vyo=(^  +  ^     -  ^)    xQ  +    (-j*  -  ^  +  ^  J  yQ    ; 

n.T=xzz=a,  ^+^ —£«.—*  +  niy  ^Brg^. 

s 

^=-tzx^o  +  Tzyy0 

,2.         32^,         32*,         321\,      32*v 

o 


"U\Z-  3x3z+  3y3z   Tz7"   3x*    3x3y  / 

L»ii  +  !!i    a2*x  +  !^x    »2».w 

"^  SySz*  3y3z   7~2"~   a_2     3z3x  /  yo 


3y      3z 

3.2  Reflection  of  waves  at  the  interface  of  two  isotropic  and 

homogeneous  media 

Let  A,  B,  and  D  represent  the  amplitudes  of  the  potential 
of  incident,  reflected,  and  transmitted  waves,  where  the  subscripts 
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%   and  s  denoting  longitudinal  and  shear  wave  (Fig.  1).   Accord- 
ing to  Huygen's  principle,  the  phase  velocities  for  each  wave 
are  equal  at  the  boundary,  therefore, 

C1A   .  Cls    .  Cl&   =  ^ls   m  C2i        m   C2s 
p  "  sine,   '  simr,    '  sine,'   '  slnff       sine2  "'  sinf2 

It  is  shown  that  the  longitudinal  or  vertically  polarized  wave, 
ty  =$     y      ,  is  always  reflected  and  transmitted  in  the  modes  of 
longitudinal  and  vertically  polarized  shear  waves  (Brekhovski Kh, 
1960) .   In  other  words,  the  vector  displacements  potentials  in- 
volved can  always  be  assumed  to  be  in  the  direction  of  the 
y-axis  and  independent  of  y.   The  potentials  in  media  #1  and  #2 
are  assumed  as  follows: 

♦l=A»*exP  [iK1A  (xsinei+zcosei)  J  +  Bj»exp  [  iKj  (xsinOi  -  zcos6i)] 
^!=A  »exp  [iK   (xsin^+zcosK^)  lyo+Bs«exp  [iK^xsint^-zcosYj)  ]yQ 

<(>2=D  exp  [iK2  (xsine2+zcos62)  ]  , 

^2=D  exp  [iK2  (xsinK'2+zcosir2)  ]  .  (3.3) 

The  boundary  conditions  at  z=0  are 

3$     3*      3$     3* 

— =-  +    *    =  — -  +   cY.  (3.4) 

3z     3x      3z     3x  |J**' 

-  -  =%-  m    -  -  =2-  (2    5) 

3X     3Z      3X     3Z  VJ,:" 
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PI6.I.    Z&Fl&C7\0l\  OP   PLKU&  WAVk  AT  rUKJB-    &OUU0AZY 
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Is  2s 

2»u5i? +  »'(i^  -  ^i=  2"^ +  »2(^  -  ^y  •  «3-7» 


Upon  substituting  Equation    (3.3)    into  Equations    (3.4),    (3.5), 
(3.6)    and    (3.7),    the  boundary  conditions   are  expressed  in  terms 
of  potentials   as   follows: 

(ArBl)Kucos61+(As+Bs)Klssinri=DJlK2jlcose2+DsK2gsiny/2    ,  (3.8) 


(Ai+Bt)Kli8ine1-(A8-B8)Kl8C06f1-DiK2l8ine1-D8K28C08jr2    J  (3.9) 

v 

pl  (A*+V  ( l-2-¥*L*2*i)  +Pi  (As~Bs)  sin2^i  " 

ls  2 

K 

P2DJl(l-2-|^sin2e2)+p2D8sin2/2      (3.10) 
K2s 

K2 
pI*A£~Bi*    -2-"sin2e1"Pi^A>+%^cos2,f'i  = 
Is 


P2D£  ^2~-sine2""p2DsCOs2^2      *  (3.11) 


K2s 

By  setting  As ,  or  A^  equal  to  zero,  corresponding  to  longitudinal 
or  shear  wave  incidence,  respectively,  the  amplitude  of  the 
generated  waves  can  be  expressed  in  terms  of  the  amplitude  of  the 
the  incident  wave  potential.   It  is  to  be  noted  that  when  the 
incident  angle  is  small,  all  the  sinusoidal  term  tends  to  the 
limits  1  and  0.   Then  the  boundary  condition  is  further  simplified 
to 

<WKirD*K2l    »  <3'12> 
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-(VBs)Kls=-DsK2s  ;  (3'13) 

-Pl(At+Bt)-p2Dt  ;  (3.14) 

"pl(VBs)=-psDs   •  (3'15) 

Solving  Equations  (3.12)f  (3.13),  (3.14),  and  (3.15) 

B£   p2CzrplCl*    D*    2p2C2£       pl 


At   P2C2£+P1CU  '  kt        p2C2A+plCl*   p2 


fs   ?2C2s-?lCls    D_s  £p2C2s     ^1 

As  "  psC2s+plCls'   As  "'  p2C2s+plCls   p2 


(3.16) 


(3.17) 


Equations  (3.16)  and  (3.17)  imply   that  when  the  angle  of 
incidence  is  small,  no  change  of  mode  occurs  at  the  interface 
of  acoustic  wave  motion.   This  is  a  very J important  limitation 
in  the  acoustic  simulation  problem.   The  angle  of  incidence  has 
to  be  very  small  if  change  of  mode  is  to  be  avoided  in  the 
experiment. 


4.   THE  GENERAL  KIRCHHOFF  SOLUTION  FOR  SCATTERING 
FROM  ROUGH  SURFACES 


4.1  General  solution  for  surface  with  one  dimensional  roughness 

Beckmann  (Beckmann,  1963)  has  derived  the  solution  for  the 
mean  scattered  field,  power,  and  the  statistical  distribution  of 
those  quantities  by  the  Kirchhoff  approximation  method.   The 
principal  limit  of  the  approximation  is  that  the  surface  must 
not  contain  a  large  amount  of  sharp  edges,  sharp  points  or  other 
irregularities  with  small  radii  of  curvature.   The  criterion  for 
the  validity  of  the  approximation  is  given  as 

4Kr  cos6>>\  , 
c        ' 

where  r  is  the  radius  of  curvature,  9  is  the  local  angle  of 
incidence,  and  A  is  the  wavelength  of  the  wave. 
The  rough  surface  is  given  by  the  function 

5  -  5(x)  (4.1) 

with  mean  level  coinciding  with  the  plane 

z  =  0  .  (4.2) 

The  medium  in  the  space  z>£(x)  is  assumed  to  be  isotropic  in 
which  a  monochromatic  plane  longitudinal  wave 

El=eikl,r-ia)t  (4.3) 


is  transmitted,  where 


K 


=il  .  Kl  (4.4) 


1  A 


K 


1 


is  the  wave  number  of  the  incident  wave,  which  is  assumed  to 
lie  in  the  xz-plane  (Fig.  2) ,  and  f   is  the  radius  vector 


r  -  xxo+5(x)zo  m  (4.5) 
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The  angle  of  incidence  is  denoted  by  B^,   the  scattering  angle 
by  62,  where  the  magnitude  of  k2  and  k1  are  equal: 

[*2I  -  I^J  =  K  -  X1  •  (4.6) 

iL  also  lies  in  the  xz-plane  for  a  one-dimensional  rough  scat- 
tering surface. 

In  order  to  deal  with  plane  scattered  waves ,  the  observa- 
tion point  P  is  removed  to  the  Fraunhofer  zone  of  diffraction, 
R'-k»,   where  R'  is  the  distance  from  P  to  a  point  B(x,£(x))  on 
the  rough  surface  (Fig.  2).   In  other  words, 

K2R'=K2Ro-K  -r  ,  (4.7) 

where  R  is  the  distance  of  P  from  the  origin.   The  scattered 
potential  E_  at  P  is  given  by  the  Helmholtz  integral 

where  + 

iKR'    iKR  -iK„ -r 
e     .  e   o   z     .  ..    Q. 

*  =  -IP-   *  R^ (4'9) 

3E 

E  and  -2—  are  the  potential  and  its  normal  derivative  on  the 
rough  surface  S.   The  values  of  those  two  quantities  are  approxi- 
mated in  the  Kirchhoff  method  by  the  value  that  would  exist  on  the 
tangent  plane  at  that  point,  i.  e., 

(E)s  =  (1+R^  ,  (4.10) 

Os  "  i(1-R)ElV"/  (4-ll} 

n  is  the  normal  to  the  surface  at  the  considered  point  B  (Fig.  3) 
and  R  is  the  longitudinal  wave  reflection  coefficient  of  a 
smooth  plane. 
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PI6.  Z  .   €>CATTSZlU(j  FZOth  A  Z01J6H   6UEFACE 


§(x; 


R6.  3.    LOCAL  eCATTEKIWS    6£0M£TCY 
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Then 


0  =  6,-8  =  9, -arc  tan?1x)  .  (4.12) 


Substituting  Equations  (4.9),  (4.10),  and  (4.11)  in  (4.8), 

(P)  =  J7)j   (l+RjE^-i  K2  nl-ifll-RlE^-n  ds 
=  J.JT  ei(V^2>'^  R(^^2)-(^2)  -S  ds, 


E. 


or 


where 


E  (P)  =  ie^)Js(E^)eiv*r^ds  (4.13) 


iL  =  K(sin6,x  -cos6,z  ) , 


iL  =  K(sin62x^+cose2z  ) , 


n     =  -sinBx  +cos8z^  . 
o      o  ' 

r  =  xx  +£  (x)  z. 


o  / 


J  =V22 


=  K(sin8,-sin60)x  -K(cos6,+cose0) z^ 

l     /  o       1     z     o 

?  ■  K(sin6,+sine0)x  +K(cos60-cos8, )  z\ 
l,     z     o       z     i  o 

ds  ■  secSdx,  ta.ft3=S'(x)  . 
For  a  one  dimensional  rough  surface  extending  from  x=-L  to  L, 
Equation  (4.13)  may  be  rewritten  in  the  scalar  form. 

«  ik^kML    (as-b)aiv  x+iVdx  (W) 

E2--nR^-j.L 
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where 

a=  (1-R)sine1+(1+R)sin62   ,        (4.15) 

b  =  (l+R)cos62-(l-R)cose1   .        (4.16) 
The  scattering  coefficient  p  is  defined  as 

p-!U    .  (4.17) 

E20 

Where  E2Q  is  the  potential  reflected  in  the  direction  of  specu- 
lar direction  (65=6,)  by  a  smooth  plane  medium-air  interface 
of  the  same  dimension.   In  this  case 


Vv=o,  5  ■  J'  ■  Oi  R  =  -1,  e1=82 


x 


so  that 

..    ikR0  f  L  ikeikRo  f  L   „ 

E20  "  ^fef—       "bdX  =  ^o~~       "2RC°SeidX    ' 
or 

ike^Lcosei        .  (4-18) 

20  1rRQ 

Hence,  from  Equations  (4.14),  (4.17),  and  (4.18), 

1 f  L(ar-b)eiVxX+iVz5dx      .  (4.19) 


4Lcos6, 

1  /-L 


For  a  smooth  surface, 


where 


p  =   1 — _  \     (1+R) cose-- (1-R) cose,   elVxXdx 
*         4Lcos6,  2  1 

J   -L 


.  _  (l+R)cose2-(l-R)cosei  .Sincv  L   ,    (4.20) 
2cos6, 


sin  VxL 
sine  V  L  ■ 


x"    V..L 
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As  a  special  case,  if  R=-l,  then  p  in  Equation  (4.20)  becomes 

p  =  sine  V  L   .  (4.21) 

o         x 

For  a  rough  surface  with  constant  reflection  coefficient  R,  a 
and  b  are  constant;  from  Equation  (4.19), 


P  = 


4Lcos 


1 


5T     "►-♦,  ,-*■-*■   T  I> 

eAV  ^dx-  —  e 
-L  z       U  - 


=  -    -1  fl  .  2R+2Rcos(91+92)  f  L  e^dx-e  (L)  . 
4Lcos91   cos91+cos92 

The  second  termacounts  for  the  edge  effect  tends  to  zero  when 
L>>\.   Thus,  ignoring  the  edge  effect, 

p  =  L.fL  e^dx  ,  (4.22) 


2L  J-L 


where 


F  =  -  Rsec  .  =   1+cob(9i+92)  (4#23) 

1     COS9,+COS92         ' 

v-r"  -  ^|-[(sin9;L-sin92)x-(cose1+cos92)  5(x)]  (4.24) 

4»2  Rough  surface  as  random  process 

The  surface  height  £ (x)  is  assumed  to  be  a  random  variable 
assuming  values  z  with  a  probability  density  w(z) ,  the  mean 
value,  denoted  by  angular  bracket  <>  ,  is 

<Ux)>  =  0     ,  (4.25) 

and  the  mean  value  of  the  integral  for  a  stationary  random 
surface  is 
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ivr,         ivvx   iv~C  , 
e    dx   =ex   e  z  dx 


-L 


-L 


■  e 


eiv*xdx 


=  X(vz) 


e  x  dx 


-L 


(4.26) 


x(vz)  ■  \  w(z)e  z  dz  is  the  characteristic  function  associ- 
ated with  the  distribution  w(z) .  From  Equation  (4.26)  the  mean 
scattering  coefficient 


<P>  -  jg  x(vz)sinc  vxL 


(4.27) 


The  variance  of  the  scattering  coefficient  D{p}  corresponding 
to  the  mean  of  the  normalized  scattered  power  is  defined  as 


<pp*>  H 


E2.2 


20 


>  =  D{p}  +  <pxp*> 


(4.28) 


The  asterisk  *  denotes  the  complex  conjugate,  and  from  Equation 
(4.22), 


-L  rh 


pp*  = 


41/ 


.iVxCxi-x^)  eivz(5l-52)  dxidx2 


'-l  ;  -l 

Denoting  eivzU1~K2)      &Y  X2tvz~vz)'  then  from  Equation  (4.28), 

L  /-L 


pp 


*      - 


41/ 


l-L 


X2CvZr-vz)    eivx(xl~x2)dx1dx2 


-L 
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"h"e  x,(v,-v  ,  .    eiv  «i-«i>    "  L  L  •t'i.-a>-iT-Ul"la>-i-a 


'2NVz  vz 
and 


C,=  S(x,)  ,   e,*?(x,)   • 


X9(v  ,-v  )  is  the  two-dimensional  characteristic  function 

£a  Z         Z 

of  the  distribution  w(  z1    ,  z2)  .  The  random  rough  surface  is 
assumed  to  have  a  Gaussian  distribution 

1       2  /0  2 

.    -z  /2o^ 


w(2)  =  /T^2     e  <4-30) 

and  the  two-dimensional  distribution  is 


w(z,,z,,t)  =  g         -exp 

1   ^      2iro^/l-c2(T) 


2  2 

z,  -2c(t)z,z2+z2 


2a2(l-c  2(t)J 


(4.31) 


2     2 

a  ss  <z  >  is  the  variance  of  the  rough  surface, 

<zlz2> 

c(t)=    *   is  tne  autocorrelation  function  of  the  rough 

<  z,  > 

surface, 

t=x,-x2  is  the  separation  factor  between  two  points  x, 

and  x~« 
The  autocorrelation  function  usually  assumes  the  following  two 

forms: 

2  .    2 

—  x  /T 

(i) .   c(x)=e   ,    gaussian  correlated,  (4.32) 

—  I  T  I  /T 

(ii) .   c(x)=e  '  "    exponentially  correlated,      (4.33) 

T  is  called  the  correlation  distance,  which  is  much  less  than 
L  to  assure  a  random  surface. 

The  one-  and  two-dimensional  characteristic  functions  are 
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evaluated  as  below  (Appendix) : 


2tto 
1    2      2 
=   e"  ¥   Vz 


1_.  e"22/20^1 


exvz"   dz 


where 


2      2 

g  =  a  v2    5 


X2(VZ'-V2)     = 


-  x    iv.,(z,-z0),      -„ 

w(z,,z2t)e     z>   1      2   dz^z2 


(4.34) 


(4.35) 


—  0L»     I  —00 


-   exp   -g(l-c(x))       .  (4.36) 

From  Equations    (4.22),    (4.26),    (4.28),    and    (4.29),    the  variance 

of   p    is 

D{p}    =    <pp*>-<pxp*> 
„    r  L  r  L 


4l/ 


eivx(x1-x2)    X2(vz,-vz)-x(vz)x*(vz)    dXldx2, 


'-L  >-L 


Using  the  relative  coordinate  x  already  defined  as 

T  =  Xx-X2   , 

and  introducing  the  center-of-mass  coordinate 

*o  =  7(X1+X2) 
the  equation  (4.37)  can  be  rewritten  as 


(4.37) 


(4.38) 


(4.39) 


D{p}=^        dx 

4l/i-L      °;-2L 


[C1J       ■  ^     3(x,,x,) 

o     lm    *  x2(t)-xx*]   jjzhj-^ 


2       r2L 

-£-.,. 2L.        [eiv*T      Xz(T)-x.X*l-det 
4L  1-2L 


1      -1/2 
1      +1/2 


dx 


2L 


=2T)_2fVxTL*2<T>-XX*>T     • 
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From  Equations  (4.34)  and  (4.36),  D{p}  can  be  expressed  in  terms 
of  c(t)  as 

D{p}=F!e-^2L  eivxT(e*cCT)-l)dT    .       (4.40) 
2L 

/  -2L 

Two  kinds  of  correlation  functions,  Gaussian  and  exponential, 

are  considered  for  g<<l?  and  g>>l  . 

(i) .   Gaussian  correlated  surface, 

(A)  If  g<l,  egc^T'can  be  expanded  by  a  uniformly 

convergent  series. 

2L 


-2   f         »  m  m,_v 

«/  i  f  --g\   «lvxt  Z.  g  c  ULa, 

2L    J-2L     m  X  m* 


If  g<<l, 


=F2  -g\   ivxx  £  gmcm(T)dT,  since  L>>T 
2L6   J _®    m=l  ml 

2    2 

_/ttTF       -g  ^  _g__  e— j .  (4.41) 

2L        e     m=i  m!/5T  4m 


D{p}  =  /l|F!g.  exp  ^-Vx2T2/4  ^   ;  (4.42) 


g*l 

_^2  e-Vvx2T2/4<    D{p}    <   £**!      .  (4.43) 

(B)    g>>l, 

,2  r2L      tmm /-,_.-^2/T?,    .1  dx 


zr2L  -t2/t2» 

.  £\        eivxTe-g(l-e       /Jtl2p2e-g 
2L   J-2L 


sinc2v  L 


|!reivxTe-g(i-e-T2^2;dt 
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Because  g  is  much  greater  than  unity,  the  integrand  is  neglible 
if  t/T  is  not  in  the  close  vicinity  of  zero;  the  integral  is 
therefore  simplified  as 


D(p%  = 


,ivx^-g[i-(i-^2/T2)]dT 


21"  2   2 

=  |_   eivxTe^T  /T  dx 


TF   ir 

2iT?exP 


2  2 


(4.44) 


4g   J 

(ii) .   Exponentially  correlated  surface:  the  same  procedure 
in  (i)  can  be  followed. 
(A)  g<l, 

.2 


D{p} 


^?*   e-g  2  SL^  /St 


_m   -„2 
ml"  "XT 


(4.45) 


m*l   m  +vx  T 


if  g<<l, 


{p }  =  — t—  •  *J— 


D{p 


L   1+vx^2 


(4.46) 


(B)  g>>l, 


F        ciT 

Dtp)  =  —  •  TTX  2  ' 

g*+vx  T 


(4.47) 


4.3  Statistical  distribution  of  the  field 

The  mean  scattered  field  and  power  is  show/?  in  Equation 
(4.2);  the  probability  distribution  of  those  quantities  can  be 
found  by  looking  at  the  random  variable  p .   From  Equation 
(3.22), 
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-L  3 


j-l 

=  re1*  .  <4-48> 

A  =  £M  is  a  constant,  and  $  .  is  a  random  quantity.   Beckmann 
L  D 

has  shown  that  for  a  very  rough  random  surface  (g>>l,<p>=o) 
with  zero  mean  <p>=o,  the  distribution  of  the  resultant  phase  \|> 
is  uniform 

wU)  =  i-   ,  -»<♦<»,  (4.49) 

and  the  distribution  of  the  amplitude  r  is  Rayleigh  distributed 


2 

P 
The  variance  of  pp*  is  found  to  be 


w(r)  =  §f-}  •"*  /°{p}  .  (4.50) 


D{pp*}=\  r4p(r)dr-D2{p} 

'0 

r-  2r5  e-r2/D<p>dr_D2{p} 
-)o  dTp"} 

.  D2{p}  (4.51) 

The  standard  deviation  of  pp* 

/D{pp*}  =  D{p}  (4.52) 

The  normalized  scattered  power  pp*  is  found  having  its  mean 
value  and  standard  deviation  equal  to  D{p} 

4.4  The  effect  of  absorption  in  the  medium 

The  presence  of  absorption  of  the  space  in  which  the  wave 
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is  transmitted  has  a  great  influence  on  the  mean  scattered 
powers.   The  damped  wave  equation  has  a  general  solution 


„    ik«r (l-ia/k)-iwt 
hi  —  e  t 

which  implies  a  change  in  the  value  of  V"x  and  V"z  of  Equation 

(4.13) . 

Let 

V  =  V  -  i^  V   , 
vx   vx   XK   vx  ' 


Vz  "  Vz  "  *t  Vz  ' 


x{vn    =.J_(   e-z2/2a2Qi(vz-iaVz)Zdz 
Z    tt) 


eAVVz  *F 


Hence, 


2   2.,     r  2   2     2 

=  e10  vz  a/k.e-|5_|z_(i-°L_j 


/   \  * i      \  ~g (l-a  A  ) 

x(v2)x*(vz)  =  e  yv    ' 


(4.53) 


(4.54) 


The  two-dimensional  characteristic  function  x„  iv",-v')  is 

X(vW)  -  ^Br.^«<W-Vk)Cl-o)-i.i-^di, 


2o' 


(Zj-CZj) 


6"  2a^(l-c2)  eivz(l-ia/k)(zr«2)d(   c  , 
00  /2tto2(1-c2)  1    2 


2«  2 


ig°A  -  |(l-o  A  \   ei(l-c)vz(l-ia/k)z2 


^2\o2 


.  eigocA>e-g(i-a2A2)(i-c)# 


,  e"(Z2/2a)dz. 


(4.55) 


-28- 


D{p} 


.  |!\   eivxT   eigac/k.exp  -g(l-a2/k2)(l-c)  -exp  [-g  (l-a2/k2)]dx 
2L  )-2L 


For  a  very  rough  surface  (g>>l)  ,  c(t)  is  significant  only  for 
very  small  value  of  x/T, 

D{P)=f  Vvx  .exp  [-gd-aV^ld-OJe^^dT  , 


or 


D{p}  =  |_     eivxx  exp[.g(i.a2/k2)(i-c)]  dT  .      (4.56) 


From  the  analogy  between  (4.56)  and  (4.40),  the  variance  of  p, 
i.  e.  D{p}  is 

(i) .   Gaussian  correlated  surface: 


Dtp}  -  -j 


5 T~  '  exP 

L   g(l-a2/k2) 


2  2 
V  'l 
x 

2 — 2~ 

4g(l-aVk' 


(4.57) 


(ii) .   Exponentially  correlated  surface: 


n/  \  -   p2   g(i-*2/k2) 
D{p}  =  L ^  "  2„2 


(4.58) 


L    gZ(l-a2A2)2TvV 


Equations  (4.17),  (4.57),  and  (4.58)  shows  that  the  abosrption 
will  reduce  the  value  of  <EE*>  ,  and  high-frequency  wave  causes 
more  attenuation  than  a  low-frequency  wave. 
4.5  Limitations  of  the  general  Kirchhoff  method  applied  to 
acoustic  wave  scattering 

The  acoustic  wave  reflection  coefficient  is  a  function  of 
the  local  angle  of  incidence,  or  alternatively,  is  a  function  of 
6,  ,  and  Ctx) . 
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R  =  R(e)  =  R^  -  arctan'U)    .      (4.12) 
The  approximations  of  <a5'>   and  <b>  in  Equation  (4.14),  which 
are  always  used  for  varying  reflection  coefficient  case,  are 
much  more  complicated  in  the  acoustic  wave  scattering  problems. 
Owing  to  the  elastic  property  of  the  medium  when  6^  increases, 
the  generated  waves  change  mode  from  one  to  another  and  the 
reflection  and  transmission  coefficients  change  very  rapidly. 
If  the  rough  surface  is  chosen  to  be  a  gently  rolling  surface 
(T>>a) ,  the  probability  distributions  of  the  slope  w(C')  are 
(i) .   Gaussian  correlated  surface: 


w(S')  -===  exP 

/4tto^ 


2,12 


4a 


(4.59) 


(ii) .   exponentially  correlated  surface: 


(4.60) 


wU')  in  Equations  (4.59)  and  (4.60)  imply  a  more  dense  distri- 
buiton  in  the  neighborhood  of  zero  slope.  If  8,  is  also  chosen 
to  be  small  angle,  the  value  of  R  can  be  averaged  over  as 

R(e1,5')  "  R(C'=<£>)=  RO^  , since  <£,'>   =  o   .  (4.61) 

For  longitudinal  wave  incidence  with  9,=0,  from  Equations  (3.17) 
and  (3.16),  the  reflection  and  transmission  coefficients  are 
defined  as 

Reflection  coefficient 

R12(61)  =  Z2"21  ,  (4.62) 

Z2+Zl 
Transmission  coefficient 


D,,(e,)  =  ;r%-&  >  (4.63) 


Z2>Z1  ?' 


I 
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where 

z      .   P±1L   ,      lf.  !l2li        .  (4.64) 

z2        cose2    '        X     cosSj^ 

If  the  #2  medium  is  free  space  or  air  and  the  #1  medium  is 
liquid  or  solid 

RlB^  =  -1   , 

RO^  =  0    .  (4.65) 

Under  the  previous  assumptions,  a  and  b  in  Equation  (4.14)  are 
considered  to  be  constant  and  the  general  Kirchhoff  approximation 
can  be  applied  to  the  acoustic  wave  problems. 


5.   ACOUSTIC  WAVE  SCATTERING  FROM  LAYER 

5.1  Layer  with  rough  interface  in  the  back 

The  evaluation  of  the  scattering  field  from  a  very  thick 
layer  with  roughness  in  the  back  is  a  direct  extension  of  the 
scattering  of  a  rough  surface.   Thickness  d  (Fig.  4)  is  assumed 
to  be 

d  >>  L, 
so  that  the  smooth  boundary  of  the  layer  is  at  the  far  field  of 
the  random  rough  boundary.   The  layer  to  be  discussed  is  a  layer 
without  absorption,  so  that  the  wave  can  transmit    in  it  with- 
out attenuation.   The  smooth  interface  has  the  following  char- 
acteristics, at  z  =  d: 

|E3|  =  D12  |El|     ,  (5.1) 

|E2I  =  D21  |E4|     ,  (5.2) 

where  E.,  E,,  E-,  and  E4  are  shown  in  Fig.  4,  and  D12,  D21  are 
the  transmission  coefficients. 
1   2Z2 


12  "  m  z2+z1 


2Z1 
D„,  «■  m 


21  "   z2+Zl 


Z2  "  P2C2^ 

Zl  =  P1C1« 

n2 


(5.3) 


m  = 

p 


1 


The  rough  surface  in  the  back  has  the  same  statistical  properties 
as  before.   Then, 


-Tl- 
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Ex  -  eikl*E+ik3dsGce3     ,  (5.4) 

E3  =  D12eiJt3-?+ik3dsec93     ,  (5.5) 

The  mean  scattered  field  <E4>  at  a  distance  R^   from  the  origin  is 

..    _    iic^r",  +ik-dsec8, 

lk-Le      4    13              3  -a/2        P 

<E„>   =  — = D-.cos  e,e  y/     •  «f  sine  v  L,    (5.6) 

4                               irR,                              12  J                   lit                  x 

where  g  and  vv  is  related  to  k-  and 

jL  =  k3  (sin  63xQ  -  cos  93  zQ)     , 

2P 

c,  sin  6,  sin  6- 

n  ^  =  lp  =  — , i.  =  — , -=-  is  the  Snell's  refraction  index.   The 

12    c,  sin  8,  sin  8. 

2p  3          4 

mean  scattered  power  for  a  very  rough  surface  (g>>l)   is 
k  2L2    ,   TF2    2  rr~  vx2T2 

<=4E4*>  =  T7  C°S  93*2lT  D12  /     9   'eXp[-  "Tg"  ]  '  (5'8) 

it  R, 


After  passing  through  the  smooth  interface,  the  mean  scattered 
field  and  power,  from  Equations  (5.2),  (5.6),  and  (5.8),  are 

ik,LeikRo+ik3d(8pce3+8pcea)  _a 

<V   "  WR+2dsec^) ^D12D21COse3'2L*e   9«^V    ' 

o  J 

(5.9) 

k    2L2D.    2D,    2COS28.    TF2       ^r—  v   V 

<E2E    *>    =   -| 12      21      2 i-S-   /-|  exp  [-   -f— ]     .(5.10) 

2   2  tt2(R  +2dsec8.)2  2L      /     g  4g 

o  J 
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Since 

ikLelkROcose1 

E2o  =         inr  ' 

o 


2    2        2 

k  L  cos  e, 

^20=20**    =  TT2 

w   Ro 


so  that,  for  R  >>  2dsec63    , 

<E~>  .  ^,.  j cose. 


<E~>  .  ~,   ,    „    COS6-   _    _ /- 

2    ,-,  rv  ^   i2k,dsec6-,     3  »F.-g/2  .__,,  T 
<P>  =  "^20  =Vl2D21e    3     3(  ScTe^  )2Te    sincvxL 


(5.11) 


Using  Equation  (4.37),  the  variance  of  p  can  be  calculated  for 
different  statistical  properties  of  the  rough  interface. 


(i)  Gaussian  correlated  surface: 

2 

~r     i  2„      2„      2rCOS    93,    TFf      /it  avrs    ,       VX    *    , 

D{p}    =   n12    D12    D21    [ 2~ ],2L-   /g  6Xp    [ 5g~]       ' 


9  2  2 

cosz6,  m„2      ,—  v„t* 


1 


(5.12) 


(ii)  .  Exponentially  correlated  surface 


2a 

cos  e,  „2 


nr    x  2n     2n     2,  "3X  P*    gT  ,-  ,,* 

D{p}  =  n12  D12  D21  ( -r).—.—^-T-  .      (5.13) 


cos  e1  g  +vx  T 


For  the  back-scattering  from  a  layer  against  air 
©2  ■  ~®i  '  ^4  =  "^3     '     R  =  -1 


i+cos(e3+e4)     2 

F  ■  -RSpc90« r — ; r =  sec  8. 

H   3  cose.+cose.         3 


v     =  k3(sine3~sine4)    =   2nilksin63   =  2ksin61  , 

vz  ■  -k3(vne3+coseJ   =  2nl2kcose3  =  -2k/n12^-siniJe1     ,(5.14) 
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2  2      2  2     2     2 

g  =  v  o  =  4k  a  (n12  -sin  6.^) 


112  "  C2* 


The  quantities  in  Equations  (5.11),  (5.12)  and  (5.13)  show  the 
backscattering  from  the  rough  surface  in  the  back  of  the  layer. 
The  backscattered  field  from  the  smooth  interface  in  the  front 
of  the  layer  is  negligible  except  for  the  normal  incidence  case, 
which  is  equivalent  to  the  reflection  in  the  specular  direction. 
If  the  pulsed  signal  is  used  to  approximate  a  monochromatic  plane 
wave,  the  backscattering  from  the  front  surface  can  be  easily  re- 
jected by  adjusting  the  gate  position  in  the  experimental  measure- 
ment. 

In  the  actual  case,  all  layers  are  more  or  less  absorptive. 
For  a  layer  constructed  with  a  material  which  has  an  attenuation 
factor  a,  <p>  and  D{p}  should  be  modified  as  follow: 

_   _    i2niokdsec0,   iy2/k,  , 3*2 

<p>  =  n12D12D21e   12      3-e     ^(ScTe^  x 

2 
|!--exp[-2dasece3-  ^(1-^-j)  ]  sincvxL     .   (5.15) 

k3 
(i) .   Gaussian  correlated  surface: 

2    2    2  cose3  2  Tp2  j — — — — 

D{p}  =  n12  D12  D21  (15551q)     -2L"  /gd-d^A^) 

2  2 
v  Tz 

•  axt>  [-4Adsec6, = 5-]     ;     (5.16) 

3      4g(l-dT/k3^) 
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(ii) .   Exponentially  correlated  surface: 

222  cose3  2  f2   gd-q2A32) 

D{P>   =n12D12D21     ^     .-.^—f-^—^. 

exp[-4adse63]  .  (5.17) 

2    2  2    2 

The  quantity  a  A3   \s  usually  very  small,  so  that  (1-a  /k3  )**1# 

but  exp[-4adsee,]  is  very  important,  it  is  nearly  unity  at  low 
frequency,  and  decreases  very  fast  as  frequency  increases. 
For  a  layer  with  a  very  rough  interface  in  the  back 

<p>  -  o    .  (5.18) 

It  is  observed  that  the  wave  backscattered  from  an  absorptive 
layer  with  rough  side  in  the  back,  the  attenuation  by  the  absorp- 
tion increases  as  frequency  increases.   Aside  from  the  attenuation 
by  the  layer,  the  wave  backscattered  from  the  rough  side  of  the 
layer  has  the  same  characteristics  as  from  a  rough  surface  of  the 
same  statistical  property,  except  for  a  changing  in  magnitude. 

5.2  Layer  with  rough  interface  in  the  front 

The  acoustic  wave  passes  through  the  rough  interface  into 
the  layer,  and  comes  back  through  it  after  being  reflected  by 
the  smooth  interface.   The  mean  scattered  field  is  negligible 
if  the  rough  side  is  very  rough,  but  the  mean  scattered  power, 
owing  to  the  complicated  phase  relationship,  is  difficult  to 
evaluate.   A  tentative  try  without  experimental  support  is  made 
here  to  look  at  some  aspects  of  the  nature  of  the  backscattering 
from  such  a  layer.   If  it  is  proved  to  be  successful,  the  same 
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method  can  be  extended  to  the  layer  with  both  sides  rough. 

As  a  first  step,  the  field  transmitted  into  the  layer  after 
the  wave  hits  the  rough  front  surface  is  found  by  considering  the 
field  at  Q(Fig.  5), 

v°>  -  t?  JK  -  *f >as    •  (5-i9) 


where 


ik-R  -ik-«r 
g      3   g        3  (5.20) 


Ro 


(E)g   =    (1+RJEj^ 


(5.21) 


|£       =  i(l-R)E.k\.r  •  <5-22> 

3n  1   1 


s 


From  the  boundary  conditions  of  the  acoustic  wave  reflection, 
for  small  incidence  angle  61   , 

k^l-RjE-j^  =  k3DE1    , 
P1(1+R)E1  =  P2DEi     ' 

(1+R)  =  — D  =  mD  (5.23) 

pl 

(1_R)  =   3D  =  Zifi.D  =  np  (5.24) 

Kl     C2p 

Substituting  Equations  (5.23),  (5.24)  into  (5.21)  and  (5.22) 

(E)   =  mDEj^     ,  (5.25) 

(||)   .  rtDE^.J     .  (5.26) 
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7 
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FI6.^.    TKAU5Ml55)Oki    THROUGH   A  KOUGH    lUTEEFACB 


-39- 


Substituting  Equations  (5.10),  (5.25),  (5.26)  into  (5.19), 


E3(Q)    "   47 


[mDE.f-iJiic    «n)    -   ^(inDE1ic1«n)  ]ds 


where 


n  =-sinBx  +  cosgz^    , 
o        o 


r  =  xxQ  +  C(x)zQ 


k\  =  k(sin9,x  -  cose1z1) 


jc,  =  nk(sin8-x„  -  cose^z  ) 


3~o 


3  o; 


w  =  kl  "  k3 


=  k(sine1  -  n  sine3)xQ  -  Moos^ 


n  cos63) z 


=  W  X   +  w  z^ 

X  o     z  o 


ds  =  secBdx    ,     tanB  =  £  (x) 


so  that, 


E3(Q) 


,  ik-R 
mke   3  o 

4lrR~ 


fL 


(ac'-bje1*  rdx 


-L 


(5.27) 


where 


a  =  D(sin6,  +  m  sin83)     , 


(5.28) 


b  ■  -D(cos6,  +  m  cose3) 
For  smooth  interface ,  C  =  C  =  0 


.  inkR 
_  .-.        inke    o 

E3(Q)  =  4Tr 
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belwxxdx 


-L 


or 


E3(Q) 


.  .  inkROT 

inKe -  D(cose,  +  m  cos9J  sinew  L 

2irR_  1  3  x 


(5.30) 


1_. 


For  the  direction   e3  =  arc   sin    (   -sine^^   ) 


E3(Q)    =   - 


.    ,     inkR_ 
inke  °l 

5TR~ 
o 


D(cose1  +  to/j_  _  1  sin2e 
n  i 


)    .  (5.31) 


E.(Q)  =  0    ,    otherwise. 


For  a  rough  interface, 


.  .  inkR_ 
t,  /«^  inke  o 
E3(Q)  = 4TR- 


w 


[-(b+a^  ) 
wz 


elw,rdx] 


J-L 


inkR, 


inke1   °D. l-nm+ (m-n) cos (e^-6i ) . 


T?R 


cose,  -  n  cos 8, 


e^dx 


-L 


or 


E3(Q) 


;inkR0  Q 


R    2L 
o 


e^dx 


-L 


(5.32) 


where , 


G  =  - 


,  „  l-mn+(m-n)cos  (6-,-e, ) 
inkD# 2   1 

2ir  *  cose,  -  n  cose3 


w«r  =  K  [(sine^n  sine3)  x-tcose^n  cose3H]     •     (5.33) 


-41- 

Equations  (5.32)  and  (5.33)  are  analogues  of  the  equations  for 
scattering  from  a  rough  surface.   The  mean  field  and  mean 
scattered  power  can  be  written  down  in  similar  form. 


einkRo  6  T 
<E3>  =  g— -^  e    sinewy 


(5.34) 


2  2 

If  it  is  a  very  rough  interface  (h  =  w  a   >>1) 


<E3>  ■  0 


and, 


(i)  Gaussian  correlated  interface^  2 


<E3E3*>  = 


TG 


ZLR  '  /  E 
o 


/  ir   *exp 


w  T 

X 


ih 


(5.35) 


(ii)  Exponentially  correlated  interface: 


<E3E3*>  = 


LR 


hT 

~2         T~~2 
hz+W  V 


(5.36) 


Equations  (5.32)  to  (5.36),  which  are  derived  under  the  assump- 
tion that  the  incidence  angle  to  be  very  small,  show  the  pro- 
perties of  the  transmission  through  a  very  rough  interface.   The 
next  step  is  to  find  the  effect  of  the  layer's  second  interface, 
thus  obtaining  the  overall  effect  of  the  layer  with  rough  side 
in  front.   This  layer  is  essentially  the  same  as  the  previous  one 
in  section  5.1,  except  it  has  been  turned  over  with  rough  side 
facing  the  incident  wave  (Fig.  6) .   The  field  E3  strikes  the 
lower  plane  interface  at  z  =  0  and  produces  a  reflected  field  E4 
in  the  specular  direction  only  (see  Fig.  6).   If  R1  is  defined 
as  the  distance  from  the  point  Oj  to  the  point  of  observation  Q 
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rcp 
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below  the  rough  interface,  then 


E3  " 


inkR,  G, 
e 1  1 

Rx   *2L 


iw,  «r, 
e   1  Id 


xl 


-L 


and  the  reflected  field  E4  is  related  to  E3  by  the  following 
boundary  condition: 

(E4)z— d  =  {R23E3}z=-d 

Under  this  boundary  condition,  the  field  E4  incident  on  the  rough 
surface  from  below  is 


i2k3dsec03      G 
(E4)z=Ux)=U(d-5)sece3'R23'2T 


e1Wl'rldxi  J-elK4,r2  (5.37) 
-L 


The  transmitted  field  E2(P)e3  caused  by  the  incidence  of 
(EJ   .  can  be  derived  from  Equation  (5.32),  here  the  subscript 
6-  denotes  that  the  field  E2(P)&  is  contributed  only  by  the  com- 
ponent of  E.  in  the  direction  of  k4.   Then, 


i2k3dsec03   ikRo     GiG2 
E2(P)e3=e2R0(d-5)sece3    *23'~^? 


ei(Vl+^2,?2)dx1dx2 


-L 


where 


(5.38) 


Gl  " 


G2  - 


inkLD12  l-mn+(m-n)cos (e«-e1) 


2.r, 


cose,  -  ncqs8. 


iKLD21.1-Si+f5"^)C08(fl2-e3) 
2tt  00383-1^  cos62 

n 


■44- 


w,    =  k.    -  k_ 


=  Msine.    -  n   sine3)xQ  -  McosOj^   -  n  cos63)z( 


x%  o  zl  o  , 


(5.39) 


w2   =  k4   -  k2 


-  k(n  sine3  -  sine1)XQ  +  k(n  cose3  -  cose^S. 
Vo  +  Wz2*o     ' 

?i  -  xiK  +  c(xi)?o 


r2  =  x2xQ  +  C(s+x-,)? 


2"*o 


It  should  be  noted  here  that  r^  and  r2  are  refered  to  origins  0^ 
and  0,,  respectively.  And  more -over,  since  the  maximum  value  of 
£(X)  is  much  less  then  the  thickness  of  the  layer,  so  that, 


Vp)e, 


ii2k3dsece3.eikR0  GjGj 


2R  dsece, 
o  3 


23       !    2" 
"      4IT 


-L 


ei(w1.r1+w2.r2)dXidX2 
-L 


(5.40) 


In  the  backscattering  direction,  62  =-91  ,  then 


G,    =  - 


inkLD.2   l-mn+fm-njcosteij-ej^) 
5fir  cos6,-n  cos93 


G,    =    - 


ikLD21    1-E+(^)COS(93!ei' 


2ir 


cos6,-n  cose. 


(5.41) 


$1  =  K(sin61-n   sine3)xQ-k  (cose^n   cose3)zQ 
$2  -  MsinBj+n   sine^x^-Mcose^n  cose3)zo 
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From  Equations  (5.40)  and  (5.41),  the  mean  value  of  E2(P)93  can 
be  found, 


j  oo   »  oo 


<E2(P)63>  =   2R  dsec( 


>i2kdsec93   ikR0     G\G2 


3  "^'7? 

r  LrL 


W(Z1/Z2ls)dz1dz2, 


-L 


-L 


,       (5.42 


where  s  in  Equation  (5.42)  is  the  separation  factor  between  ^ 


and  £2  ' 


s  =  (2d+C1+52)tan63 


since  iw<<d 


s  *  2.dtane. 


(5.43 


Hence,  from  Equation  (5.42) 


<E~(P)6->   =  A-exp[-h(l+c(s)]  'sinctkLfsinej^-n  sin63)] 


•sinc[kL(sin6.+n  sin63)]       (5.44) 


where 


,,2  2 
h  =  W„  a 
z 


W  =  -K(cose1~n  cos63) 


A  = 


ik^dsece-i   ikR_ 
e  3 -j  »e   o 

2dR^sec6, 
o     3 


(5.45) 


From  Equation  (5.44),  it  is  seen  that  if  the  front  interface  is 
very  rough  (h  >>  1) , 
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<E2(P) 93>  ■  0 


and, 


<E2(P)>  = 


tt/2 

< 

-ir/2 


<e2(p) e3>de3  =  0 


(5.46) 


The  evaluation  of  <E2E2*>   is  very  much  involved,  approxi- 
mations have  to  be  used  throughout  the  derivation.   From  equation 
(5.40), 


L  L  L  L 


(E2E2*)e3ss  ***" 


■*       .-*■ 


iwl*  (rl*rl)+iW2'  (r2"r2)dx1dx2dx'1dx2   (5.47) 


-L-L-L-L 


where  the  subscripts  1  and  2  are  referred  to  the  origions  0^   and 
02,  respectively  (see  Fig.  6). 

The  mean  value  of  (E2E2*) e3  over  the  rough  interface  is 


L  L  L  L 


<(E2E2*)63>  -  AA* 


<e 


iWz(51-C1,)+iWz(51-C2,>> 


-L-L-L-L 


.eiwxltxrXl)+iWx2.(x2"  x2)dx1dx2dx1'dx2«    (5.48) 


and  from  Mood  (  Mood  1963  ) ,  the  4th  order  characteristic  func- 


tion associated  with  W2  is 
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(s,T>S<eiW*(5r«i)+iWzU2-«2l> 


00   /  00  />    00   /■  00 


—  OO   —00   —00  /  —  00 


W(z1,Z2,z£,z-,T,s)ei^Ul-«i'+iW^2-«2) 


dz.dz2dz£dz2 


=  exp[-  |(2+2c(s)-2c(t1)-c(s+t1)-c(s-t1) 

-  |(2+2c(s)-2c(t2)-c(s+t2)-c(s-tz)3   ,       (5.49) 


where 


t,=x1-x^  ,  and  t2=x2-x2 


(5.50) 


For  a  Gaussian  correlated  random  surface, 

2   2 

,  ,     -t7T 
c(t)  =  e       , 

and  as  seen  from  Equation  (5.49),  the  4th-order  characteristic 

function  X4(s,t;l,t2)  is  equal  to  zero  except  when 


s<<T,  t,<<T,  and  t2<<T 
In  the  case  shown  in  Equation  (5.51), 


(s^^t^^e-^dx2-^2)/^]. 


(5.51) 


(5.52) 


Introducing  new  coordinates, 

yl  =  l^l^V  '    and  y2  =  2  (x2+x2*  ' 
and  using  Equation  (5.52)  for  a  very  rough  layer,  Equation  (5.48) 

becomes 

2  2 


(E2E*)93    =AA*« 


.LrLrL,L  t,  t2 

\\\\        e"2hV    +    ?>.«i^lT+i,teT*T1dTadyi 

4W  'exp  L~5K— J 


-L-L-L-L 
2 


4L2AA..  /£ 


Wx\t21 
exP  I  — SH- 


or 


2 

<(E2E*)9>=2L2AA**^-  .exp 
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,2  ,„2 


8h    l  ■ 


(5.53) 


where 


2  2  2 

h=k  a  (cos8.-ncos83)   , 


(5.54) 


W2  +W2  =  2K2(sin29,+n2sin  9.J 

Substituting  Equation  (5.54)  into  Equation  (5.53), 

/  ,-r,   r,*x«  \       2irT2L2AA* 

<(.E2E*)63>   =   -5-5 5 

3       K  0  (cos81-ncos63) 

•  ex;p  [-T  (sin261+n2sin2e3)/4o2(cose1-ncos63)  J 

(5.55) 
In  Equation  (5.55),  ^(E^jjpe^  is  significant  only  when  91 
and  93  are  very  small  angles,  i.  e. , 

sine, =6,   ,   cos9,=l 

sine3  =  83   ,   cos93=:l   .  (5.56) 

Then,  the  limit  in  the  integration  for  obtaining  <E2Ei>  can 
be  extended  to  infinity,  so  that 

^E2E2>=\  <<E2E2)93>d63 

«\      2L2AA*--^ y   .exp  f-T2(62+e2)/4a2(l-n)2   de, 

J..  K<r(l-n)z  L      J  J      J 

(i-mn^m-n)^  (1-^.(1-1)) 


2  2    4    4    2      2 

AA*  R23  nKLP12D21 


4d"R*»  J6L"  16ir  (1-n) 
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Rewriting  Equation  (5.55), 

2   2   2.   2       4 

<E2E*>=  R23D12D21L !*_  #k2#I  exp^ejT2/4a2(l-n)2]   (5.57) 

256*   2d2R2m2|l-n| 

2   2 

Since  E2QE*=^i2-  ,    and<E2>=o,    the  variance  of   the   scattering 

*   Ro 
coefficient  p   becomes 

<E2E2> 
D{p}  =^  *       - 

E20E20 

4 
R__D, 0D0, (1+m  )       -,  2  2    2      2  . 

=   23  12  21 x  —V-  .exp  -e£r/4a  (l-np         (5.58) 

256  Trm2|l-n|        ad 

From  Equations  (5.46),  (5.55)  and  (5.58),  it  is  seen  that  the 
backscatter  from  a  layer  with  rough  side  in  front  is  much  smaller 
than  that  from  a  layer  with  rough  side  in  back.   Moreover,  the 
parameters  of  the  media  has  more  influence  on  the  backscattered 
power  as  compared  with  the  previous  one  in  Section  5.1,  and  as 
the  incidence  angle  e,  increases,  the  backscattered  power  drop 
off  very  rapidly. 
5.3  Discussion  of  the  derivation 

(A)  Rough  side  in  back:  The  fields  E3  and  E4  in  the  layer  have 
been  assumed  to  be  caused  by  a  plane  wave  and  the  amplitudes  are 
independent  of  the  thickness  "d"  of  the  layer.  If  the  thickness 
d  is  much  larger  then  the  illuminated  length  2L,  the  amplitudes 
of  E^  and  E.  will  depend  on  the  value  of  d.  From  Equation  (5.3), 
introducing  the  dispersion  caused  by  the  thickness  of  the  layer, 
the  value  of  E3  given  by  Equation  (5.5)  will  be  changed  to 

in  keik3-ri+ik3dsece3 

E,  m — — s— LD,  -(cose,+mcose,)  .     (5.5a) 

3  ZifR-.  -L*      J-       J 
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A  similar  modification  should  be  made  on  E4  and  E2  also.   The 
results  of  these  changes  in  E2,  E3,  and  E4  will  modify  D{/>}  of 

Equations  (5-12)  and  (5-13)  to  the  form 

2   2 
fi99  0  4432r—  V     T 

n12D12D21      008^63(0086^1000863)       _kLTF     h         exp    _   x 

D(p}= t~2 2  74"    Jg  4g~ 

32*  m  cos^  d    ^y 


(5.12a) 


(ii)  exponentially  correlated  surface: 


D{p}=n12D12D21  cos263(cQseL+mcose3)4>  k4LV  _    gT      (5.13a) 

le/m2      "         cos2e1  d4  g2+vV     . 

(£)   Rough  side  in  front:   Under  the  same  condition,  D{p} 
in  Equation  (5.58)  should  also  be  modified  to  the  form 


D{o}=n2K2R23D22D21(l+m  )  .TL2.eXp|'-g2T2/4o2(l-n)2  ].    (5.58a) 

_,    3    2,,       !  74"  L  J 


l2K2R23D|^l+m4: 

T      5 

256*  m    |l-n|  od 

D{p}  in  Equations  (5.12a),  (5.13a),  and  (5.58a)  are  derived 
under  the  condition  that  the  layer  is  very  thick.   However,  in 
most  cases,  d  is  not  so  much  larger  than  the  illuminated  length, 
and  the  Equations  (5.12),  (5.13),  and  (5.58)  give  us  a  good 
estimate  on  the  value  of  D{p}. 


6.   EXPERIMENTAL  WORK 

6.1  Experimental  set-up  and  procedure 

The  experimental  measurement  of  the  variance  of  the  scatter- 
ing coefficient  has  been  done  in  the  Underwater  Acoustic  Labora- 
tory, Electrical  Engineering  Department,  Kansas  State  University. 
Details  of  the  equipment  can  be  found  in  the  report  by  Toliver 
(Toliver  1965) .   The  block  diagram  of  the  experimental  set-up  is 
shown  in  Fig.  7. 

In  the  experiment,  pulsed  signals  generated  by  the  pulsed 
oscillator  are  sent  out  and  collected  by  one  pair  of  transducers 
in  the  water  tank.   For  each  single  pulse  sent  out,  the  received 
signal  will  contain  a  train  of  pulses.   A  gating  circuit,  which 
is  synchronized  by  the  delayed  trigger  output  from  the  pulsed 
oscillator,  is  used  to  select  the  portion  of  the  pulse  train  for 
feeding  into  the  detector  and  boxcar  circuit.   The  boxcar  cir- 
cuit has  the  function  of  converting  the  discrete  pulse  into  an 
analog  signal  so  that  it  can  be  recorded  by  the  graphic  level 
recorder. 

The  distance  S  from  the  transducers  to  the  front  side  of  the 
target  is  determined  by 

S  ■  Dcose,  (6.1) 

where  6.  ■  0°,  5°,  and  10°,  and  D  is  the  distance  travelled  by 
the  radiated  signal  before  it  hits  the  target.   It  is  desired 
to  have  D  as  large  as  possible  so  that  the  illuminated  area  will 
be  much  larger  than  the  correlation  distance  of  the  random 
interface.   The  choice  of  32  inches  is  made  to  give  an  average 
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value  of  L  =  0.75  inches  for  different  pairs  of  transducers.  The 
pulse  recurrence  frequency,  PRF,  is  also  determined  by  D.   It 
is  selected  to  be  the  largest  value  without  causing  any  overlapp- 
ing of  the  first  few  returns  with  the  following  pulse  sent  out 
from  the  transmitter.   And  the  pulse  length,  which  should  be  as 
larger  as  possible  to  simulate  a  monochromatic  wave,  is  limited 
by  the  thickness  of  the  layer.   The  pulse  length  must  be  less 
than  the  time  that  is  necessary  for  the  wave  to  make  a  round  trip 
in  the  layer.   A  safety  choice  of  PRF  and  pulse  length  for  the 
model  target  constructed  and  the  distance  D  specified  is  as 
follows: 

pulse  length  =  20  *isec, 
PRF  =  250  pps. 
The  operating  frequencies  are  chosen  to  be 

0.72,  1.0,  1.28,  1.6,  1.9,  2.25,  3.0,  and  3.5  mc. 

Before  the  measurements  of  D{p}  start,  the  folliwng  pro- 
cedures are  conducted  at  8,  -  0°,  S  =  D  =  32  inches,  f  -  1  mc. 

(a)  The  target  suspension  is  carefully  checked  by  the 
returned  pulse  position  on  CRO  to  make  sure  that  D 
is  equal  to  32  inches  for  all  possible  positions  of 
transducers. 

(b)  Transducers  focusing  is  done  by  adjusting  the  mount- 
ing of  transducers  for  maximum  return  from  a  smooth 
plane  target  32  inches  apart. 

After  checking  on  the  mounting  of  target  and  transducers, 
measurements  proceed  as  follows: 

(c)  Setting  the  operating  frequency  with  the  help  of  test 
oscillator  and  CRO. 
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(d)  The  pulse  length  and  PRF  is  set  at  pulse  length  ■  20psec, 
PRF  =  250  pps. 

(e)  The  RF  output  level  and  the  transducer  compensator  is 
adjusted  to  yield  the  best  possible  undistorted  pulse. 

Procedures  (c) ,  (d) ,  and  (e)  are  conducted  iteratively  to 
fit  all  the  figures  required. 

(f)  Measurement  |e2Q|  :  the  water-air  interface  is  used 
as  target,  because  it  acts  as  perfect  reflector;  and 
D  =  32  inches,  6.,  =  0°  are  carefully  checked.   Then 
the  magnitude  recorded  by  the  recorder  gives 

M1   =  K1AB|E2()|  (6.2) 

where  A  is  the  magnitude  of  the  output  pulse,  B  is  the 
gain  of  the  transducers,  and  K1  is  the  overall  gain  of 
receiver  and  recorder. 

(g)  Varying  6,  to  the  desired  angle  (0  ,  5  ,  10  )  and 
setting  S  according  to  the  relation  shown  in  Equation 
(6.1),  then  I EJ  ,  which  is  a  function  of  the  horizontal 
position  of  the  transducers,  is  measured  by  scanning 
the  transducers  through  the  target.   The  gate  is  ad- 
justed in  the  way  that  only  the  return  from  the  back 
side  of  the  layer  is  detected  and  recorded.   The  sample 
magnitude  of  a  point  on  the  recording  sheet  gives 

M2  =  K2AB  |E2|  ,  (6.3) 

where  K2  is  the  overall  gain  in  this  measurement.   The 
ratio  of  K2  to  K.  can  be  read  from  the  settings. 
Since  <p>=0,  so  that  the  variance  of  the  scattering  coef- 
ficient is 
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.(P>-<|P|2>.<J^>   «<».(£) 


E20  Ml 


or 

.2 


n/    .        <M2>        .n    (kdb) 
D{P}  =  — Y~  *   10   WlT  / 

Ml 

where  K ,,  =  (K~/K, )  in  db  is  the  ratio  of  the  overall  receiver 

(Id      I      1 

gain  between  the  measurement  of  \E-\  and  | E2Q  \ .   Attention  must 
be  paid  not  to  overdrive  any  stage  in  the  steps  (f )  and  (g) . 
K~  is  adjusted  to  give  the  maximum  possible  recording  without 
causing  saturation  in  any  one  of  the  amplifiers  and  in  the 
recording  devices. 
6.2  Description  of  layer  target 

The  target  is  a  block  of  plexiglas  with  a  one-dimensional 
roughness  on  one  side.   (Fig.  8) .   The  length  of  the  plexiglas 
is  limited  by  the  dimension  of  the  water  tank;  the  width  is  equal 
to  6  inches,  which  is  much  greater  then  the  illuminated  area  of 
any  pairs  of  transducers;  and  the  thickness,  which  should  be  as 
thick  as  possible,  is  limtied  to  2  inches  by  the  material  avail- 
able at  the  time  of  construction.   The  profile  of  the  rough  sur- 
face on  the  layer  is  desired  to  have  a  Gaussian  probability 
distribution  in  height  with  standard  deviation  o  =  0.05  inches, 
and  to  have  a  gaussian  autocorrelation  with  respect  to  the 
horizontal  position,  where  the  correlation  distance  is  T  =  0.15 
inches.   The  profile  is  calculated  by  trial  and  error  with  the 
help  of  an  IBM  1620  computer.   It  is  noted  that  the  heights  be- 
tween -3a  and  +3o  were  considered  in  the  calculation,  and  the 
sample  heights  obtained  from  the  truncated  normal  distribution 
are  so  arranged  that  there  is  no  sudden  or  periodical  variations. 
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The  profile  obtained  was  then  cut  by  a  shaping  machine.   A  sample 
measurement  at  intervals  of  0.05  inches  was  then  made  on  the 
finished  target  and  the  distribution  and  correlation  functions 
were  calculated. 

The  measured  standard  deviation  and  correlation  distance  are 
0.0465  inches  and  0.15  inches,  respectively.   Moreover,  the 
measured  correlation  function  lies  between  Gaussian  and  exponent- 
ial, and  consists  of  certain  periodical  variation  when  the 
separation  factor  is  greater  than  the  correlation  distance  T. 
Moreover,  the  distribution  of  slope  is  found  to  be  different 
from  normal  and  to  have  a  pair  of  extra  peaks  at  £'=tan(10  ). 
This  might  be  caused  by  the  shape  of  the  cutter  of  the  shaping 
machine. 
6.3  Measurement  and  results 

The  quantity  M,  of  Equation  (6.1),  owing  to  the  drift  in 
the  electronic  circuit  and  the  disturbance  in  water  tank,  is 
usually  a  function  of  time.   After  the  disturbance  in  the  water 

has  died  away,  the  fluctuation  of  M,  with  time  is  nearly  zero. 

2 

The  quantity  M2  is  a  varying  positive  quantity,  the  mean  of  M2  is 

<  M2  >"  K  gj"2i   '  (6-5) 

and  M2i  denotes  the  ith  sample  of  M2«   The  variance  of  pp*  is 

as   follows 

2      2 

D{pp*}    =    <p    p*    >    -    <pp*> 

<M2*>    .    10(Kdb/5)-<^10Kdb/5    (6'6) 


Ml 


where 


<"a>"iJl"Si-  <6-7) 
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The  samples  {M„.}  are  taken  at  equal  intervals  from  the  recording 

sheet.   Since  part  of  target  is  not  seen  at  bigger  6^  the  number 
of  samples  taken  is  equal  to  96,  90,  and  84,  for  ex  =  0  ,  5  ,  and 
10°  respectively.   D{p}  calculated  from  the  experimental  data  is 
shown  in  Fig.  9.   As  a  comparison,  the  value  D{p}  for  backscat- 
tering  from  a  rough  surface  is  also  obtained  and  shown  in  Fig.  10. 
6.4   Discussion  of  results 

To  compare  the  experimental  results  with  the  theoretical 
solution  of  the  backscattering  from  layer,  D{p}  is  calculated  for 
6,  =  0°,  5°,  and  10°,  from  Equations  (5.12)  and  (5.13). 

The  parameters  of  water,  air  and  plexiglas  are 

(A).   Water: 

p   =1.0  gm/cm*, 

4- 
c  =  15  x  10  cm/sec, 
w 

4 
z  =  p  c  =  15  x  10  gm/cm  -sec  . 
www 


(B) .   Air: 


p   =1.29*10  gm/cm, 

a 

4. 
c  =  3.4  x  10  cm/sec, 

4       2 

z   =  0.0041  x  10  gm/cm  -sec. 
a  3 


(C) .   Plexiglas: 

p   =1.2  gm/cm, 

c  =  27,8  x  10  cm/sec 
P 

4         2. 

z  =  33.4  x  10  gm/cm  -sec. 
P 

The  reflection  and  transmission  coefficients  at  the  plane 
boundary  of  two  isotropic  medium  as  given  by  Equation  (5.4)  are 
calculated  for  the  following  cases 
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E—  x  -£  =  1.15, 


12   z  +z      p 
p  w     ^p 

2z      p 
D21  =  y   »   x  -E.  0.744 
zp+zw     pw 

z  -z 

R22  "  z  +z*        *  _1   * 
a  p 

Other  parameters  are  given  as  follows 

CW     n     ca 

n,  ~  =  —  =  0.54  . 
12   cp 

The  effective  illuminated  area  represented  by  L  is  a  rather 

complicated  quantity,  it  is  assumed  here  to  be  equal  to  the 

illuminated  area  by  one  pair  of  beam-limited  transducers;  mean 

value  is  taken  for  five  pairs  of  different  transducers. 

L  =  0.75  inch. 

The  theoretical  solutions  Equation  (5.12),  and  (5.13)  are 
calculated  and  shown  in  Fig.  11(a),  1Kb),  11(c).   The  experi- 
mental results  are  also  shown  in  these  figures. 

The  relations  among  frequency,  wavelength,  and  ratio  a/ A 
related  to  plexiglas  are  calculated  and  tabulated  below: 

Frequency  (mc)  A (mm)  Ratio  a/A(tf»1.18  mm) 

0.72  3.86  0.306 

1.00  2.78  0.425 

1.28  2.17  0.544 

1.60  1.74  0.679 

1.90  1.465  0.806 

2.25  1.235  0.955 

3.00  0.927  1.27 

3.5  0.794  1.49 

From  Fig.  11a,  lib,  lie,  it  is  seen  that  the  measurements  of 

D{p}  deviate  considerably  from  the  theoretical  solutions  from  a 

layer  having  either  gaussian  or  exponential  correlated  rough  inter- 
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face.   This  can  be  explained  by  the  existence  of  an  attenuation 
factor  a  in  the  layer.   If  a  is  a  non-decreasing  function  of 
frequency,  the  Equations  (5.12)  and  (5.13)  should  be  modified  to 
the  form  of  Equations  (5.16),  and  (5.17).   The  modified  results 
are  such  that  D{p}  is  reduced  at  low  frequencies  and  decreases 
faster  as  frequency  increases.   For  such  a  choice  of  a,  the 
experimental  results  are  then  expected  to  lie  between  the  modi- 
fied theoretical  solutions.   Since  the  target  constructed  is  such 

~r2/T2 
that  the  correlation  of  the  thickness  lies  between  e   '        and 

-  I  T  I  /T 

e  '  "  ,  the  experimental  work  proves  that  Equations  (5.16)  and 
(5.17)  are  a  good  prediction  of  the  backscattering  of  an  acoustic 
wave  from  a  layer  with  the  rough  side  in  the  back. 
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7.   CONCLUSIONS 


The  backscattering  from  a  layer  depends  on  the  thickness, 
acoustic  impedance,  and  the  statistical  parameters  of  the  rough 
side.   It  also  depends  on  the  angle  of  incidence  and  operating 
frequency.   The  thickness  of  the  layer  enables  the  separation 
of  the  second  pulse  from  the  first  one.   In  case  of  an  absorp- 
tive layer,  the  wave  propagated  through  the  layer  is  attenuated. 
The  attenuation  may  strengthen  the  frequency  dependence  of  D{p}, 
if  the  attenuation  factor  a  is  a  function  of  frequency. 

D{p}  obtained  from  a  layer  with  rough  surface  in  the  back 
has  substantially  the  same  form  as  that  froa  a  rough  surface, 
on  the  condition  that  the  layer  is  constructed  with  non-absorp- 
tive material.   If  the  random  side  is  very  rough,  the  value 
of  D{p}  decreases  very  fast  as  the  ratio  o/X  increases  where  o 
is  the  standard  deviation  of  the  rough  interface;  and  for  small 
angle  of  incidence,  D{p}  increases  as  o/X  increases. 

If  the  layer  is  turned  over  with  rough  side  facing  the 
incident  wave,  the  evaluation  of  D{p}  involves  an  integral  which 
contains  a  four-dimensional  characteristic  function  associated 
with  w  .   Approximation  is  made  by  assuming  e3  very  small  in 
evaluating  the  integral.   The  theoretical  solution  for  this 
model  has  a  very  strong  dependence  on  the  thickness  d,  even  if 
the  layer  is  non-absorptive. 
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APPENDIX:  EVALUATION  OF  INTEGRALS 


The  integration  of  the  improper  integrals 


e-x2/2a2eibxdX/ 


O      *5 

-x  /2a  cosbxdx,  and 


rx2/2a2ei(b-ic)xdx 


involves  contour  integration.   Choose  the  contour  as  shown 


2    2 

e-z'/2a'dz=lim  ^  e-x  /2a  dx+lin, 


a2b 


*-X 


e-(R+iy)2/2a2dy 


R+00 


-R 
-R 


R+OO 


R-»-oo       j  R^oo 

R 


a2b 


=    0      . 


or, 


<~ 


2/o    2  [  2/0    2 

(fee-2  /2a  dz  =\   e"x  /2a  dx  - 


-6c+ia  b)/2a    - 
e  '        dx 


s  o, 
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or 


■(x+ia2b)/2a2+  e-(x-ia2b)/2a2dx  >  \       e"x  /2a     dx 


e(^2j       e-^/2a2  e-ibx+eibx     dj£  =   f^2      . 
'  o 


-x2/2a2  w  ^  /IT      ^-a2b2/2j 

»*/'*»     cosbxdx  =  ^|  j  aev  ' 


2^2 

(1) 


Since 


-x  ?2a2   ibxdx  =    e-x  /2a  cosbxdx  +\  e"X  /2a  sinbxdx 


e  A  ''"   .e     "dx  =  e 


f°°        2         2 
=   2]    e"x  /2a  cosbxdx+O    , 
/o 


So  that 


rx2/2a2>eibxdx  =   /2ira2      .e-a2b2/2      #  (2) 


And  also 


00  2 

rx2/2a2>ei(b-ic)xdx     \        e-(^2+cx)  >eibxdx 


2      2       ..  ,         2    .         2    2 
,x-a  c*         ib(x-a  c)      a_c 
-( — * — )    .e  .e— 5 — 

e      2a2  2 

.eia2bcd(x-a2c) 

=   /27aT.eia2bc.exp  [-^(b2-c2)]    (3) 
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ABSTRACT 

The  backscatter  from  a  rough  surface  is  usually  calculated 
by  using  Kirchhoff 's  approximation.   In  this  report,  the  Kirch- 
hoff's  approximation  is  extended  to  the  backscatter'ing  of  an 
acoustic  wave  from  a  rough  layer.   The  random  rough  interface  of 
the  layer  is  assumed  to  have   one-dimensional  Gaussian  distri- 
buted surface  heights. 

Gaussian  and  exponential  autocorrelation  functions  are  used 
to  represent  the  correlation  of  the  heights  at  two  different 
points.   Expressions  for  the  variance  of  the  scattering  coefficient 
are  derived  in  the  case  that  the  rough  side  is  very  rough. 

Experimental  investigations  were  conducted  at  ultrasonic 
frequencies  on  a  target  designed  to  have  Gaussian  distribution 
and  correlation.   The  measured  variance  of  the  scattering  coef- 
ficient p,  D{p},  has  the  following  properties: 

(A)  D{p}  is  highly  frequency  dependent;  it  decreases  as 
the  frequency  increases. 

(B)  D{p}  decreases  more  rapidly  than  that  obtained  from  a 
rough  surface. 

(C)  The  dependence  of  D{p}  on  the  incidence  angle  8,  is 
such  that  it  increases  as  9,  increases,  if  6,  is  small. 

The  experimental  results  lie  between  the  theoretical  results 
calculated  for  two  kinds  of  correlation.   This  agrees  with  the 
measured  value  of  the  correlation  function  of  the  target. 


